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NATURAL NONEQUILIBRIUM STATES IN QUANTUM STATISTICAL MECHANICS. 

by David Ruelle*. 



Abstract. A quantum spin system is discussed, where a heat 

O , flow between inhnite reservoirs takes place in a hnite region. A 

Q ! time dependent force may also be acting. Our analysis is based 

CN I on a simple technical assumption concerning the time evolution 

^ ' of infinite quantum spin systems. This assumption, physically 

natural but currently proved for few specific systems only, says 
that quantum information diffuses in space-time in such a way 
that the time integral of the commutator of local observables 
^ I converges: J dt \\[B , a''^ A]\\ < oo. In this setup one can define 

i> ■ a natural nonequilibrium state. In the time independent case, 

this nonequilibrium state retains some of the analyticity which 
O ■ characterizes KMS equilibrium states. A linear response for- 

^ ! mula is also obtained which remains true far from equilibrium. 

Q\ [ The formalism presented here does not cover situations where 

^ ■ (for time independent forces) the time translation invariance and 

,-Q ■ uniqueness of the natural nonequilibrium state are broken. 

-(— > 
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Introduction. 

Traditional nonequilibrium statistical mechanics has been focussed on approach to 
equilibrium (Boltzmann and followers) and on situations close to equilibrium (Onsager 
reciprocity, Green-Kubo formula). More recently, a fruitful rigorous study of nonequi- 
librium steady states for classical systems far from equilibrium has been initiated, using 
the concept of Gaussian thermostat [6], [11]. Among the results are the Gallavotti-Cohen 
fluctuation theorem [8], [9], the Dettmann-Morriss pairing rule [2], [16], and a general lin- 
ear response formula [13] (see also Dorfman [3], Ruelle [14] for reviews). In the approach 
just referred to, flnite classical systems are driven away from equilibrium by nonhamil- 
tonian forces, and cooled by a Gaussian thermostat. The more natural approach which 
uses Hamiltonian forces and inflnite heat baths is more difficult, and results there are still 
preliminary [4], [5]. 

Compared with the classical theory, quantum statistical mechanics exhibits significant 
differences: equilibrium states (KMS states) are more intrinsically tied to the dynamics, 
and the forces are fundamentally Hamiltonian. In particular, the use of a Gaussian ther- 
mostat does not appear feasible. We are thus led to studying inflnite systems with Hamil- 
tonian forces. Fortunately, the dynamics of inflnite quantum spin systems is relatively 
amenable to study: a C*-algebra A is associated with the system, and the time evolution 
is described by a one-parameter family (a*) of automorphisms of A. 

The physical situation which we wish to discuss is that of a flnite quantum system E 
interacting with inflnite reservoirs Ra, {a — 1,2, . . .) in equilibrium at different tempera- 
tures, chemical potentials, .... The system S is also acted upon by a force which may be 
time dependent: 



Time 

dependent 

force 



For notational purposes it is convenient to write T, = Rq. Note that a flnite system 
subjected to time dependent forces will in general heat up, and that a single reservoir Ri 
can act as a thermostat. Another case of interest is when two reservoirs Ri, R2 at different 
temperatures interact via E (no time dependent force is assumed here). Variations of the 
setup just described have been considered by a number of authors (see in particular Hepp 
and Lieb [10], Spohn and Lebowitz [15], Jaksic and Fillet [12]); an important feature of 
the present approach is that it uses more realistic reservoirs. 

The strategy of analysis that we shall adopt in this paper will be to compare the 
time evolution (a*) of the interacting system described above with a noninteracting time 
evolution (a*) where E, Ri, R2, . . . evolve independently. Let a be an invariant state for 
(a*), where the reservoirs i?i, i?2, • • • are at temperatures /3i , /3^ , . . .. It will be possible 




to determine a nonequilibrium natural state pt for the interacting system by the condition 
that it reduces to a in the distant past. 

To estabhsh the desired relation between the evolutions (a*) and (a*), we assume that 
the interactions between E and i?i, i?2, • • •, and the time dependent force acting (possibly) 
on E are of local nature. There is however at this point a serious technical problem: our 
definition of natural nonequilibrium states requires that, when A, B are spatially localized, 
time integrals of the type j dt \\[B,a*A\\\ converge (and similarly for integrals involving 
(a*)). These conditions (further discussed in Section 2.1) are physically natural but have 
been proved to hold only for very special quantum spin systems. 

In view of the above difficulty, we shall in what follows adopt an axiomatic approach. 
We shall make technical assumptions (A1)-(A5) on the dynamics of our quantum spin 
system, and derive our results from these assumptions. While (A1)-(A4) could easily 
be seen to hold for specific systems (see [1] Section 6.2), there is a problem with (A5), 
as discussed above. The interest of the results obtained seems however to justify our 
axiomatic approach. 

In Section 1 we discuss (A1)-(A5) and derive the existence of a *-isomorphism ujt 
between the C*-algebra A of the full interacting system, and the C*-algebra Ay of the 
union of the infinite reservoirs i?i, R2, ■ ■ ■ (E being omitted). The isomorphism ut inter- 
twines between the time evolution (a*) of the full iteracting system and the noninteracting 
time evolution (a*) restricted to the union of the reservoirs (see Theorem 1.6). 

In Section 2 we show how to define nonequilibrium natural states pt asymptotic in the 
distant past to noninteracting equilibrium states {pt is related to these states via uJt). In 
Section 3 we consider the case of time independent forces, and assume that the equilibrium 
states of the noninteracting reservoirs are KMS states at different temperatures. The 
nonequilibrium state pt = p is now time independent and retains some of the analyticity 
which characterizes KMS states. 

In Section 4 we obtain a linear response formula for dp when there is a small change 
dh in the interaction of the finite system E with the infinite reservoirs (the temperature 
of the reservoirs is not changed). This quantum linear response formula holds far from 
equilibrium and is very similar to the corresponding formula for classical systems [13]. 
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1 Assumptions (A). 

The following assumptions (A1)-(A5) specify our mathematical setup. 

(Al) A, and Aq, Ai, A2, ■ ■ ■ are (Gnitely many) C*-algebras with unit elements such 
that A is the C* tensor product of the Aa (a = 0, 1, 2, . . .), and Aq is the algebra ofnxn 
matrices for some finite n > 0. 

There may be several norms on the tensor product ^a>oAa satisfying the C* property 

||A*A|| = ||A|p and the cross-norm property || (S>a>o ^a|| = na>o ll^all- -'-^ ^^^ -^a are 



identified to operator algebras on Hilbert spaces Ti^, the operator norm of (S>a>0'^a acting 
on ®a>o^a does not depend on tfie cfioice of tlie faitliful representations, see [1] Section 
2.7.2. Tfie closure ^ayo-^a of (S>a>0'^a for this operator norm is our C* tensor product of 
the Aa- We shall later also use Ay = ^a>oAa and denote by Iq, 1> the unit element of 
A, A>. 

The algebra ^o corresponds to the system T, = Rq, and Ai, A2, . . .to the reservoirs 
Ri, R2, . . . , while A describes the total system. 

(A2) (a*), (a*) are one-parameter families of *-automorphisms of A such that 

a* = ®a>0Oia 

and (al) is a one-parameter group of *-automorphisms of Aa, i-e., a^J = identity, «„«„ = 
al+\ iaiAar=aiA:. 

We shall write a(t, s) = a*~*, a(t, s) = {a^)~^a^, and also a>(t) = a> = (S>a>o «!• 

The evolution (a*) describes the uncoupled systems T^ = Rq, Ri, R2, ■ ■ ■ , while (a*) 
describes the total system where E is coupled to the reservoirs and subjected to time- 
dependent forces. 

(A3) There is a dense subset V C A and for each s G R there is e > such that, if 
D eT), the functions 

t' ^ a{t',t)a{t,s)D , t' h^a{t\t)a{t,s)D 

are (norm-)differentiable when t,t' G (s — e,s + e). 

Clearly, one can assume that "D is a *-subalgebra of ^. Note also that the families (a*), 
(a*) are strongly continuous because for each s G R and A E A, the functions t 1-^ a^A, 
a^A, or equivalently t 1— *> a~^A, a~^A, or equivalently for each B E A the functions 
t 1-^ a{tj s)B, a{tj s)B are continuous on (s — e, s + e) as uniform limits of continuous (in 
fact differentiable) functions 1 1-^ a{t — s)D, a{t, s)D with D eT). 

(A4) There is a finite dimensional linear space T such that 

A ® i> c J^ c ^ 

and a function h : H ^' J-' such that h is bounded continuous, self-adjoint (h = h*), and 
with the notation of (A3) 

-^ia{t', t)D' - a{t', t)D') |,,=, = -i[h{t), D'] 

if 

tG(s-e,s + e) , D' e a{t,s)VU a{t,s)V 

This expresses that the interaction between the system E and the reservoirs is of local 
nature. 



1.1 Lemma 

If A E A, the functions 

t h^a^a^y^A, a^a~^A 

are differentiable with derivatives 

t ^ -ia^[h{t), (a*)"^A] , ia^[h{t), a~^A] 



We first take A = aW, D eV. Then (A3) impUes that the function t i-^ a*{a^)~^A 
is differentiable for t G (s — e, s + e) because 

(a*'(a*')-^^-«*("*)~^^) =«*'^ (a(t',t) - a{t' ,t))a{t, s)D 

Using (A4) we see that the derivative is 

-ia*[/i(t), a(t, s)D] = ~ia^[h{t), (a*)"^A] 

so that 

a\aY^A = a%a')-^A - i f rfwa"[/i(w), («")"^A] (1) 

J s 

and this formula remains true for all A G ^ and all s, t G R [cut [s, t] into small subintervals 
and use the density of V in ^]. Since u — > a'^[h{u), {a'^)~^A] is continuous, we see that 
the derivative of t i-^ a^{a^)~^A is t i-^ —ia^[h{t), {a^)~^A]. The case of t i-^ a^a~^A is 
similar, with 

a^a-^A = a'a-'A + i f dua''[h{u),a-''A] (2) Q 

J s 

Note also the formulae 

d ^ 

—a{s, t)a{t, s)A = —ia{s, t)[h{t), a{t, s)A] (3) 

—a{s,t)a{t,s)A = ia{s,t)[h{t),a{t, s)A] (4) 

which follow directly from Lemma 1.1. 

1.2 Lemma 

Let a^ be obtained from a by the replacement h -^ h + Xk. Then the functions 
X -^ Q;^(t, s)A are differentiable and 



d 
dX 



J s 



Writing A A = A' — A and Act = 0;'^ — a'*', we obtain from (4) 
— Aa(s, t)a(t, s)A/AX = iAa{s, t)[/i(t), a(t, s)A]/AX + ia^' [s, t)[k{t),a{t, s)A] 

(Jib 

and using again (4) this is 

= Aa(s, t)a^(t, s)/AA— a^(s, t)a{t, s)A + ia^' (s, t) [A;(t), a(t, s) A] 
Writing $(s,t) = q;^(s, t)a(t, s), $^(s,t) = Aq;(s, t)a(t, s)/AA we have thus 

(-$^(s,t))$(s,t)-i$(s,t)A=$:^(s,t)$(s,t)-i-$(s,t)A + za^'(s,t)[A;(t),a(t,s)A] 

hence 

(^$:,(., t))$(., t)-i$(., t)A + $:,(., t){^Ms, t)-^ms, t)A = za^' (., t)[fc(t), a(t, .)A] 
at at 



hence 

($^(s,t))$(s,t)"^)$(s,t)A = m^ (s,t)[/c(t),a(t, s)A] 



dt 

hence 

(a^' (s, t)«^(t, s)A - A) /AX = {^'^{s, t))$(s, t)-^A 



i / ctra^'(s, r)[/c(r),a^(r, s)A] 

J s 



(5) 



From this it readily foUows that $/\^(s, t)A has a hmit $'(s, t)A when AA -^ 0, so that the 

j_ 



derivative -jT-a^{s,t)A exists (and is equal to $'(s, t)a(s, t)A). Thus 



d /■* 

(— -a^(s,t))a^(t, s)A = z / (ira^(s,r)[A;(r),a^(r, s)A] 

dA y^ 

or 

— -a^(s,t)A = z / dra^{s,T)[k{r),a^{r,t)A] 
dX J,, 

as announced. [] 

Our last assumption will play a crucial role. 

(A5) There are dense subsets S G A and £y C Ay such that, if F G J-', E G S, 
Ey G £y, then 

I ds\\[F,{aT^E]\\<^ 

J —QO 

ds\\[F,a-%lo®Ey)]\\<oo 



00 




We may take for £ (resp. £y) a *-subalgebra of ^ (resp. Ay). Condition (A5) may 
be taken to mean that, as s ^ — oo, a^E, q;*(1o®-E'>) diffuse rapidly away from a bounded 
region, in particular the region of interaction between the system E and the reservoirs (see 
the further discussion in section 2.1). 



1.3 Proposition. 



There are *-niorphisms u~[ : A ^ A and oo^ : Ay -^ A such that, for ah A E A, 
A> eA>, 

lim a{t,s)a{s,t)A = ujf A 

s — y — oo 

lim a{t, s)a{s,t) {Iq^ Ay) = uj^ Ay 



Using (1) and (2) we see that 

a'{a')-^A^{a^)-^A + i / dta'[h{t),{aY^A] 

J s 

a^a~%lo0Ay) =a°(lo® A>) -i / dta^[h{t),a~\lo Ay)] 

J s 

Since h{(j) is bounded and takes values in the finite dimensional space JF, (A5) shows that 
the right-hand sides converge when s — > — oo provided A — E & £, Ay = Ey G £y. By 
density of £^ in .4. and £y in Ay , the limits 

uj^A= lim a''{a'')~'^A , Uq Ay = lim a'^a'^'ilo ® Ay) 

s — > — oo s — > — oo 

exist for all A E A, Ay E Ay, defining *-morphisms Uq : A ^ A and Uq : Ay -^ A. 
Therefore the limits asserted in the Proposition also hold, with u^ = a~^UQa^, u^ = 



or 



1.4 Proposition. 

If A E A, uj^A = lo ujfA where ut is a *-morphisni A -^ Ay. 
Notice first that if E E £, (A5) gives 

/ ds sup \\[X^ly,{a')-^E]\\<oo 
/ ds sup \\[Y®ly,a'{a')-^E]\\<oo 

J-oo YeAo,\\Y\\<l 



The derivative of s i-^ a''^{a^) ^E has bounded norm in view of Lemma 1.1, and the 
function 

s^ sup \\[Y®l>,a\a')-^E]\\ 

YeAo,\\Y\\<l 

has thus bounded Lipschitz constant. Therefore 

hm sup \\[Y®l>,a\a')-^E]\\=Q 

or 

(VreA) lira [Y®l>,a'{a')-^E] = Q 

s — > — oo 

Putting a^ {a^)~^ E in matrix form (eij(s)) with z, j G 1, . . . , n, we see that eij{s) ^ if 
i ^ i and eu^s) — ejj{s) -^ when s -^ — oo. Writing 



1 "^ 



n 

i=l 



we have hms^_oo \\a''^ {a^)~^ E — Iq ® -E'(s)|| = 0. For general A, we can also write 
a^{a^)~^A in matrix form {aij{s) and define A(s) = - Yl7=i '^"(^) ■ Approximating A 
hy E ^ S shows that 

lim lla'^'aM- lo® A(s)|| = 

s — » — oo 

or equivalently that uj^A = Iq ® UtA where Ut is a *-morphism A -^ Ay. Q 
1.5 Proposition. 

Ut : A ^' Ay and u^ : Ay -^ A are reciprocal * -isomorphisms. 

Choose t G R and A E A. Writing B = UfA, A' = iv^B, we can in view of Propositions 
1.3, 1.4 choose 5" < such that if s < 5", 

\\a{t, s)a{s,t)A - lo (^ B\\ < e (6) 

\\a{t,s)a{s,t){lo®B)-A'\\<e (7) 



(6) implies 

and (7), (8) give 



\A-a{t,s)a{s,t){lo0B)\\<e (8) 



\\A-A'\\ <2e 
hence u^UtA = A' = A. If Ay G Ay, a similar argument gives UtUJ^ Ay = Ay. [] 
We may summarize our findings as follows 



1.6 Theorem. 

There is a *-isoniorphisni uJt : A ^ Ay such that 

lim a{t, s)a{s,t) A = 1q ® iVt A 

s — > — oo 

lim a{t, s) a{s, t) {1q ® Ay) = uj^^ Ay 

s — > — oo 

In particular uJta{t, r) = a>(t, T)uJr- 

This follows from Propositions 1.3 and 1.5. [] 

2 Nonequilibrium states. 

We shall call natural nonequilibrium states those states which, for the evolution (a*), 
reduce in the distant past to equilibrium states for the evolution (a*). This definition is 
possible because in the distant past (a*) and (a*) are close to each other as a result of our 
assumption (A5). We now discuss further this assumption. 

2.1 Discussion of assumption (A5). 

Physics suggests that, when A, B are spatially localized, the estimate 

||[S,a*A]||^t-^/2 

typically holds for a (i-dimensional quantum lattice system. Specifically, some examples 
treated by Bratteli and Robinson ([1] 5.4.9 and 6.2.14) conform to this diffusive type of 
behavior (these examples are however rather special). When (i > 3, the t"'^'^ estimate 
implies that 



j dt\\[B,a^A]\\ <oo 



which is basically our assumption (A5). Consider now the case of time independent forces, 
i.e., let (a*) be a one-parameter group of automorphisms. If we have 



/oo 
dt\\[B,a*A]\\<oo forA, Bg^ 
-oo 



where £^ is a norm-dense *-subalgebra of A, Bratteli and Robinson say that (a*) is L^{S) 
asymptotically ahelian ([1] Definition 5.4.8). Under this condition they prove the existence 
of our *-morphism uj^ ([1] Proposition 5.4.10) which they call M0ller morphism by anal- 
ogy with quantum scattering theory. Bratteli and Robinson point out that the difiiculty 
in proving L^ asymptotic abelianness in particular models is not surprising since the exis- 
tence of the M0ller morphism is a form of ergodicity*. The approach of [1] has technical 



* I am indebted to Derek Robinson for pointing out Section 5.4 of [1] in connection 
with our assumption (A5). For an example of nontrivial study of ergodicity in an infinite 
system, the reader is referred to Fidaleo and Liverani [7]. 
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advantages over the approach adopted here in Section 1, but our discussion has the interest 
of applying to time dependent forces (and of being self contained). 

As we have said, the assumption (A5) means that a^E or a''^{l (E) Ey) rapidly diffuse 
away from the region of interaction between the system E and the reservoirs Ri, R2, ■ ■ ■ 
Such a diffusion is possible because the reservoirs are infinite, and more precisely of di- 
mension > 3. This dimensional restriction is physically not surprising if we think of a 
macroscopic description of the state of our system by a continuous temperature function T 
tending to finite values (3^ , /9^ , ... at infinity in the different reservoirs. In the simplest 
case T should satisfy the heat equation AT = 0, but if rf = 1, or 2 this implies that T is 
constant or unbounded. We are thus forced to imagine that our reservoirs have dimension 
3 or more. 

Physically one expects that (for time independent forces) the time translation invari- 
ance and uniqueness of the natural nonequilibrium state may be broken. The formal- 
ism presented here does not cover these situations. Indeed our natural nonequilibrium 
state p will, under natural assumptions*, satisfy the mixing property limt^oo \p{Aa*B) — 
p{A)p{B) I = and therefore p has no nontrivial decomposition into time invariant or peri- 
odic states. In conclusion, we expect that some but not all situations far from equilibrium 
are covered by our assumption (A5). 

2.2 Definition of natural nonequilibrium states. 

Let a a be a state on Aa, invariant under the one-parameter group (a*) for a = 
0, 1, 2, . . . We shall later impose that the a^ with a > satisfy the KMS condition (see 
below). The GNS construction gives for each a a Hilbert space Ha, a representation it a of 
Aa by operators on Ti^, a vector O^ such that 

aa(-) = (Oa,7ra(-)^a) 

and one-parameter groups Ua{-) of unitary operators such that 

Ua{t)na = 0„ , Ua{t)na{A)Ua{t)-^ = naia^A) 



* The state p will turn out to be conjugate to the product ®a>oO"a of KMS states 
describing reservoirs, and (using central decomposition) it is natural to assume that the a^ 
are factor states. Therefore ®a>oC"a and p are factor states. Assume also the asymptotic 
abelianness condition 

lim 1 1 [A, a*S] 1 1 = for A, B e A 

t — >oo 

(this is implied by L^ asymptotic abelianness). Then the mixing property 

lim \p{Aa^B) - p{A)p{B)\ = 

t— >oo 

holds (see [1] Example 4.3.24). 

10 



In particular a = ®a>oC"a is a a^-invariant state on A. We say that the time dependent 
state pt on ^ is a natural nonequilibrium state (NNES) if it is of the form 

Pt = lim a{s,t)*a 

s — > — oo 

or 

Pt{A) = lim a{a{s,t)A) = lim a{a{tj s)a{Sjt)A) 

s — > — oo s — > — oo 

= a{uj+A) = ay{uJtA) 
where a> = ®a>oC'"a- We may thus write 

Pt = WiV> 

which shows that the NNES pt does not depend on the initial state ctq of the system E. 
Our definition gives in particular 

Pt{a{t,T)A) = pr{A) 

We also have 

pt{A)=a{A) + i dua{[h{u),a{u,t)A]) (9) 



[where we have used the formula — z J^ dua{t, u)[h{u), a{u,t)A] = A—a{t, s)a{s, t)A, which 
follows from (3)]. 

2.3 The KMS condition. 

Let Pa > 0. The (a* ) -invariant state aa satisfies the /S^-KMS condition if, whenever 
A, B E Aa, there is a bounded continuous function F on {z : < Imz < /3a}, analytic for 
z < luiz < Pa and such that for all real t 

aa{B.a'A)=F{t) , aa{a'A.B) = F{t + iPa) 

[The physical meaning of this condition is that a^ is an equilibrium state at temperature 

Pa']- 

We shall say that pt is a /?-NNES if it is a NNES associated with /3a-KMS states aa- 
It describes thus nonequilibrium in the presence of reservoirs Ra at various temperatures 

Pa'- 

We assume for simplicity that the Wa are faithful representations [this is natural: it is 
physically reasonable to assume that the Aa have quasi-local structure, with simple local 
algebras, so that the Aa are simple algebras. See [1] Section 2.6.3]. 

3 Time independent forces and nonequilibrium steady states. 

We shall now consider the situation where the forces acting on the system E are time 
independent: h{t) = h, a(s,t) = a*~^ , Ut = UJ, and the NNES pt = p is a nonequilibrium 
steady state (NESS). We have thus 

UJO*' = CXyU , p = u*ay 

11 



p{A)=a{A)+i f dua{[h,a-''A]) 

J — oo 



p{B.a^A) = a>{uB.ai,ujA) (10) 

when A, B E A. If the o"a are /3a-KMS states for a > 0, we shaU say that p is a /3-NESS. 
Consider the elements B G Ay such that 

||S||i = nminf{VTT||Sj-„|| : ||S - V ®„>oSj-a|| < e} 

j a>0 j 

These elements form a *-algebra A^ with norm ||.||i > ||.|| 

Let 7i> = (8>a>o^a and i3(7i>) be the *-algebra of bounded operators on 7i>. Since 
the 7Ta are faithful, the map ^a>oT^a '■ ^a>oAa -^ BiTiy) extends to a faithful *-represen- 
tation 7r> of Ay by bounded operators on 7i>. 

Write t = (ti,t2,...) and let d* = ®a>Qofa = «(^i) ® 0(^2) ® . . ., i.e., d* is the 
automorphism of Ay such that 

ny{a'A) = (H t/.(t„))7r>(A)(n Uy{ta))-^ 
a>0 a>0 

(d* is unique because 7r> is faithful). 

3.1 Proposition. 

If A E Ay, B G A^, there is a complex function F of the complex variables Za, 
continuous and bounded by ||A||.||i?||i on Ylayoi^a • ^ — ^^^^a < f^^a}, analytic in 
Ylayoi^a '■ < Im^a < (3a} , and such that 

ay{B. a^ A) = F{t) , ay {a^ A.B) =F{t + i/3) 

where P = (/3i, /32, . . .). 

Take first A, B in the algebraic tensor product ®a>o-^a, viz., 

i j 

Then 

ayiB.a'A) = J2Il ^a(S,ad*-A,„) 

ij a>0 

extends, by the KMS condition for the aa, to a function F bounded and continuous on 
na>o{-2^« • - ^^^a < Pa}, and analytic in na>o{-^« • < Im^^ < Pa}- Using the Cauchy 
formula in several variables we also see that |F(2i, Z2, ■ ■ ■)\ has the sup-norm 

||F|| = max snp\F{ti + i'r]iPi,t2 + ir]2P2,---)\ 

VI, V2,--- t 

12 



where ?7i, ?72,- • • take the values or 1. Using the KMS condition and separating the indices 
a' with I3a' = from the indices a" with fia" = 1 we find that 

F(tl + Z?7l/3l, t2 + Z?72/32, • • •) = XI W'W(^a'{Bja'-0^'^A^a')(ra"{a^'-A^a'.Bja") 



^ ^a{®a'Bja'{®aCi^'"Aia) ®a" Bja") = "^ a {(g) a' B j a' A^ A (g) a" Bja") 
ij J 

hence 

iiFii<ii^iiEnii^^-»ii^ii^ii-ii^iii 

Using the density of ®a>oAa in Ay and in Al, concludes the proof of the proposition. [] 

Note that Proposition 3.1 remains true if one changes the assumptions to S G A^, 
AeAl- 

3.2 Corollary. 

There is a dense *-subalgebra A^^^ of A such that if A e A, B e A^^^ or B e A, 
A E A^'^\ the function t i-^ p{B.a^A) extends to a continuous function on {z : < Imz < 
m.ma(3a}, analytic in {z : < Imz < m.ma(3a}- 

Taking A^^^ = uj~^Ay, this follows from (10) and Proposition 3.1. [] 

Under suitable physically reasonable conditions one should be able to take A^^' — 
A^ = {A e A : \\A\\i < oo} where 

||A||i = liminf{VTT||A,„|| : ||A - V ®„>oA,g| < e} 

J a>0 j 



3.3 The modular group of p. 

As pointed out by M. Winnink*, if (Tip, tt^, O^) is the cyclic representation associated 
with p, then O^ is cyclic and separating for TVp{A)", and therefore a modular group (r*) 
of automorphisms of TTp^A)" is defined. In fact we may write Tip = 7i> = ®a>o'^ai 

TXp = TTyUJ = (®a>07ra)w, Op = 0> = (ga>0^a and 

Therefore the modular group (rt) corresponds asymptotically in each reservoir Ra to the 
noninteracting evolution (a^) accelerated by the factor —(3a- 



* Private communication. 
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4 A general linear response formula. 

4.1 Proposition. 

For a perturbation dh{-) of the time dependent interaction h{-), the time dependent 
nonequihbrium state p. satisGes the following linear response formula 

Spt{A)=i drpt{[a{t,T)6h{T),A]) 

J —QO 

when A G {a^)~^S. More precisely, if /c : R ^ JF is bounded continuous self-adjoint and p^ 
is the NNES corresponding to the interaction h{-) + Xk{-), then A ^^ Pt{A) is differentiate 
at X = when A G {a^)~^S, and 

d /■* 

—p^{A)\x=o = i J drpt{[a{t, r)k{r),A]) 

Using (9) we have 

p^{A)=a{A) + i lim / dua{[h{u) + Xk{u),a^{u,t)A]) (11) 

where a^ is obtained from a by the replacement h \-^ h + Xk. From (5) we get also 

a^{u,t)A-a^{u,t)A = i / rfr a^(w, r)[A/c(r), a°(r, t)A] (12) 



where a^ = a. From (11) we obtain 

p^{A)-p',{A)= lim (Ai(.) + A2(.)) 



where 

ft 



Ai{s)=i dua{[Xk{u),a^{u,t)A]) 

J s 

A2{s)=i / dua{[h{u) + Xk{u),a^{u,t)A-a'^{u,t)A]) 

J s 

and, using (12), (3), 

A2(s) = - / du I dr(T([/i(w) +A/c(w),a^(w,r)[A/c(r),a°(r,t)A]]) 

J s J u 

/"* r d 

^ -i I dr I du—a{a^{u,r)[Xk{T),a'^{T,t)A]) 



s J s 
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= -z / rfr((T([AA;(r),aO(r,t)yl])-a(a^(s,r)[AA;(r),a°(r,t)A])) 
so that 



(p,^(A)-pO(A))/A = z lim / dTa{a\s,r)[k{r),a\r,t)A]) (13) 

If A G (a*)~^£^ we can, given e > 0, choose sq such that 

rfT||[/c(r),a(r,t)A]|| < e 



if s < Sq. Therefore (13) imphes that A t-^ Pti^) is differentiable at 0, with 

d /■* 

p^{A)\x=o = i lim / dra{a{s,r)[k{r),a{r,t)A]) (14) 

s — > — oo / 

J — oo 



rfA 



= i hm / (ir(j(Q;(s, t)[Q;(t, r)A;(r), A]) 

s — > — oo / 

J — oo 

and the proposition foUows readily. [] 

4.2 Corollary. 

Under the conditions of the Proposition 4.1, for time independent h{-) and p. = p, we 

have 

d /■* 

—p^{A)\x=o = tJ dTp{[a{t,T)k{T),A]) 

^oo 

= i / dspila'^'kit- s),A]) 

Jo 

4.3 Remarks. 

If one can choose £ independent of A such that 

dT\\[kiT),a\T,t)A]\\ 

D 

converges when A E £, uniformly for A G (Ai,A2), then A -^ Pti^) is continuously 
differentiable on (Ai,A2), and 



^Pt{A)=if drp^,{[a\t,T)k{T),A]) 



Higher order derivatives are also defined, and 

d"^ 
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p^^([fc(ri),a^(ri,r2)[Mr2),---[Mrn),«^r„,t)A] •■■]]) 
as is seen by repeated differentiation of 

d /■* 

p^(A)=i lim / (ircr(a^(s,r)[A;(r),a^(r,t)A]) 
dX s^-ooJ_^ 

(obtained by restoring the dependence on A in (14)), and Lemma 1.2. We may also write 

^" p^(A) = %^n\ [ ■■■ [ dTn---dTi 



pU[a\t, ri)A;(ri), [a^(t, r2)A;(r2), ■ ■ ■ [a\t, r„), A] ■ ■ ■]]) 
For time independent /i(-) and k{-) = k, we have thus 

1 Jn /"OO POO 

--— p^(A)|A=o = ^" / dai--- dar,p\[Ka\-ai)[kr ■ ■[k,a\-ar,)A]- ■ ■]]) 
ni ciA Jq Jq 
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